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Abstract 

The thermal expansion and contraction of ring particles orbiting a planet or 
asteroid can cause secular orbit evolution. This effect, called here the thennal expansion 
effect, depends on ring particles entering and exiting the shadow of the body they orbit. A 
particle cools off in the shadow and heats up again in the sunshine, suffering thermal 
contraction and expansion. The changing cross-section it presents to solar radiation 
pressure plus time lags due to thermal inertia lead to a net along-track force. The effect 
causes outward drift for rocky particles. For the equatorial orbits considered here, the 
thermal expansion effect is larger than Poynting-Robertson drag in the inner solar system 
for particles in the size range -0.001 - 0.02 m. This leads to a net increase in the 
semimajor axis from the two opposing effects at rates ranging from -0.1 R per million 
years for Mars to — 1 A per million years for Mercury, for distances -2 R from the body, 
where R is the body’s radius. Asteroid 243 Ida has -10 R per million years, while a 
hypothetical Near-Earth Asteroid (NEA) can have faster rates of -0.5 R per thousand 
years, due chiefly to its small radius compared to the planets. The thermal expansion 
effect weakens greatly at Jupiter and is overwhelmed by Poynting-Robertson for icy 
particles orbiting Saturn. Meteoroids in eccentric orbits about the Sun also suffer the 
thermal expansion effect, but with only ~0.0003e 2 AU change in semimajor axis over a 
million years for a 2 m meteoroid orbiting between Mercury and Earth. 
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1. Introduction 

Radiation forces (photon thrust) can cause the secular orbital evolution of small 
particles orbiting planets and asteroids. One of the best-known such forces is the 
Poynting-Robertson effect, which causes orbit decay for particles circling the Sun 
(Poynting, 1903; Robertson, 1937; Bums et. al, 1979). Similarly, it also causes orbit 
decay for ring particles orbiting a planet. 

The Yarkovsky effects (diurnal, seasonal, and Yarkovsky-Schach) can also result 
in the secular orbital evolution of ring particles; for instance, those which compose 
Saturn’s rings. A previous paper (Rubincam, 2006) dealt with how photon thrust from 
Yarkovsky thermal forces may affect the orbital evolution of Saturn’s main rings (see 
also Vokrouhlicky et al., 2007). While it was found that Yarkovsky orbital evolution 
could be significant over the age of the solar system assuming no destruction of particles, 
the conclusion was tempered by the fact that all of the various Yarkovsky effects depend 
on the linkage between insolation (from the Sun or Saturn) and the orbit: particles must 
be in undisturbed principal axis rotation for the Yarkovsky mechanism to work optimally. 
Any disturbance in the linkage, such as from collisions, weakens Yarkovsky. Particles 
tumbling on timescales much shorter than the orbital period would tend to average out the 
effects to the point where Y arkovsky forces would become insignificant. 

The diurnal and seasonal Yarkovsky effects also vanish at the opposite extreme of 
synchronous rotation; here there is too much linkage between orbit and insolation. Hence 
the Yarkovsky forces operate on a restricted set of rotational states. 

The present paper examines yet another photon thrust effect, called here the 
thermal expansion effect. Like the Yarkovsky-Schach effect, it depends on particles 
entering and exiting the planet’s shadow for its operation. And like the Yarkovsky- 
Schach effect, the thermal expansion effect tends to drive stony particles away from the 
planet, helping keep rings alive, thus working against forces like Poynting-Robertson and 
the seasonal Yarkovsky effect, which make ring particle orbits decay. Unlike the 
Yarkovsky effects, it is largely blind to spin state insofar as a particle’s radial temperature 
distribution is concerned: particles can be in principal axis rotation, tumbling, or in 
synchronous rotation and still suffer the thennal expansion effect. Hence the effect is not 
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hampered by Yarkovsky’s restrictive spin conditions. The effect does not work well in 
densely populated rings; neither does Yarkovsky. Therefore the thermal expansion effect 
will operate best in thinly populated rings, such as those around Jupiter and conjectured 
at Mars (e.g., Hamilton, 1996). 

The basic idea behind the thennal expansion effect is qualitatively illustrated in 
Fig. 1. A spherical ring particle revolves in a circular orbit about a planet. Solar radiation 
pressure pushes on the particle, as denoted by the thick arrow. The direction of the force 
is away from the Sun, and the magnitude of the force is indicated by the size of the arrow. 
The force is proportional to the particle’s cross-sectional area, which for a sphere is jtR 2 , 
where R is the particle’s radius. Just before entering the shadow, the solar radiation 
pressure pushes on the particle in the direction of motion, which increases the semimajor 
axis. When the particle dives into the shadow the solar heating turns off, and the particle 
cools; thennal contraction causes the radius R to shrink. When the particle emerges from 
the shadow, solar radiation pressure acts against the motion so as to decrease the 
semimajor axis. However, the radius R upon emergence is smaller than when it went into 
the shadow, so that the solar radiation pressure is less. The Sun does warm the particle 
back up, increasing the radius, but this takes some time due to thermal inertia. As a result, 
the forces do not balance when averaging around the orbit, and the net effect is a positive 
along-track acceleration which increases the particle’s orbital semimajor axis. The 
particle moves away from the planet. 

The effect depends on both thennal expansion and contraction for its operation; it 
is called here the thermal expansion effect for short. It depends crucially on the particle’s 
thennal inertia. If the particle heats up and cools down instantaneously and thus expands 
and contracts instantaneously, then the radii R at various points around the orbit is 
symmetrical with respect to the shadow, and the effect vanishes. At the other extreme, a 
particle with large thermal inertia will not heat up or cool down much, leading to little 
expansion or contraction and a small effect. Hence somewhere between the extremes the 
effect works best. 

The quantitative treatment of the thermal expansion effect follows. It is based on 
the formalisms of Rubincam (1987, 1998, 2006). The ring particles are in circular orbits 
about the planet. The particles are assumed to be made of rock in the case of Jupiter and 
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planets and asteroids sunward of Jupiter, while for Saturn the particles are composed of 
water ice. In all cases the particles are taken to be opaque spheres with homogeneous 
properties. The particles are assumed to be far enough apart from each other so that the 
shadowing or illumination of one particle by another can be ignored; in other words, the 
rings are sparsely populated. 

The present paper concentrates on rocky equatorial rings around the terrestrial 
planets (though none have been observed) and Jupiter (which definitely has rings), and 
the water ice rings of Saturn. It also looks at two asteroids: 243 Ida, and a hypothetical 
NEA (Near-Earth Asteroid). Uranus, Neptune, and Pluto are ignored, due to the weakness 
of the effect thanks to their great distance from the Sun. The ring particle orbits are 
assumed to remain circular; changes in orbital eccentricity e are not considered. Readers 
uninterested in the mathematics can skip to section 6. 

2. Insolation on ring particles 

The ring geometry is shown in Fig. 2. The Sun is in the x-z plane. A particle orbits 
in the planet’s equatorial plane in a circular orbit at a distance a from the planet’s center. 
Its angular position (f) in the equatorial plane is measured from the x-axis. The insolation 
incident on the particle can be written as a Fourier series 


is the semimajor axis of the Earth’s orbit, and api an is the semimajor axis of the planet’s 
orbit about the Sun. For Mars, for instance, ap/ an = 1.52 AU. Also, (j) = nt, where t is time 
and 



2 


( 1 ) 


where Fp = 1361 W nf 2 is the solar flux at the Earth (Kopp and Lean, 2011), ctp = I AU 



( 2 ) 
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is the mean motion of the particle, with Mp/ an being the planet’s mass and G being the 
universal constant of gravitation. 

The coefficients bj express the fact that sunlight turns off when the particle enters 
the planet’s shadow and turns back on when it exits the shadow. It will be assumed here 
that the Sun is a point source, and that there is no refraction of sunlight by the planet’s 
atmosphere into the dark region behind it. In other words, the transition from full sunlight 
to complete darkness is instantaneous; thus (1) expresses a “boxcar” function. In this case 

=1-/1, (3) 

and 

bj = 2sin [/'(l- A)jr\/jt, (4) 

where A is the fraction of the orbit which is shadowed (e.g., Rubincam, 2006), so that A 
ranges from 0 (orbits outside the shadow) to 1/2 (orbits hugging the planet). 

3. Ring particle temperatures 

The variable insolation will cause the temperature T of the particle to fluctuate. 
The temperature distribution inside the particle can be found from the heat diffusion 
equation 


KA 2 T = pC p — 

p dt 

(e.g., Chapman, 1967), where K is thermal conductivity, p is density, and C p is specific 
heat; K, p, and C p are taken to be constant throughout the particle. Temperature T can be 
written T= To + AT, where To is a constant in space and time. Only the radial temperature 
distribution inside the particle will be considered here, so that AT = AT(r), where r is the 
distance from the center of the particle. The above equation then reduces to 
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K 



r 1 dr 



( 5 ) 


Switching to complex fonn, the time-dependence of Fs un can be written in terms of 
exp(z'vf), where i = (-1) and vis the frequency of the /th term in (1), so that v=jn. It is 
well-known that for each tenn in the insolation (1) there will be a corresponding term in 
the temperature of the form ATj = constant x jo(kr) exp (ivt), where 


j 0 (kr ) = 


sin kr 
kr 


( 6 ) 


is the spherical Bessel function of order zero, and 


k = 



If the particle radiates like a grey body, then the thermal emission will vary like T 4 . The 
constant in AT can be found from linearizing the power balance at the surface r = Ro, 
where Ro is the radius of the particle, using the approximation T 4 ~ To + 47o 3 AT, so that 


K 


d(AT) 
dt , 


+ 4e IR oT 0 3 AT = (l-A v )F Sun /4 


Here Ay is the particle’s albedo in the visible, £m is the infrared emissivity, and a is the 
Stefan-Boltzmann constant. The factor of 4 on the right side comes from the solar flux 
being incident on jtRo and being averaged over the particle’s surface area of AjtRo. 
Writing AT = 2 / ATf r), by the above equation 
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ATj(r) = - 


(1 -A v )h.F £ R 0 

f a E 

zjp(kr) 

4 K 

\ a Plan 

) ^zcosz + c 0 sin 


with z = kRo = x(l - i ) and 


2 K 


,1/2 


with c 0 being 


K 


-1 . 


At the surface of the particle a typical temperature term by (7) is 


A7VR 0 ) = 


(1 A v )bjF E R 0 ( a E j B t +B 2 i(vt-Sj) 


4 K [a 


Plan , 


B 


where 

zj 0 (z) B l -iB 2 
zcosz + c 0 sinz fi 3 

with 

B\= x (sin 2x + sinh 2x) + Co (cosh 2x - cos 2x), 
Bj = x (sinh 2x - sin 2x), 


(V) 


( 8 ) 


( 9 ) 


( 10 ) 


( 11 ) 

( 12 ) 


Bt, = 2x 2 (cosh 2x + cos 2x) + 2cox (sinh 2x + sin 2x) + cq 2 (cosh 2x - cos 2x), 


(13) 


Rubincam 


4/25/13 


10 


and 


tan <5 ; = — . 
7 B x 


( 14 ) 


4. Ring particle thermal expansion and contraction 

Because of temperature changes inside the particle, the molecules at r will be 
displaced radially a distance u{r) due to thermal expansion and contraction. Elasticity 
theory says u{r ) is given by the expression 


u(r) = 


a v (l + cr p ) 

' 1 

9(1 -a P ) 

r 2 


J r 1 AT {r) dr 


+ 


2r(l-2cr p ) 
^(1 + CT P ) 


f\ 2 AT(r)dr , 


(15) 


where ay is the volume expansion coefficient and Op is Poisson’s ratio (Landau and 
Lifshitz, 1970, p. 22). The radius of the particle will be given by R = Ro + AR. At the 
surface (15) becomes 


A * = “W-^/VA7-(,-)rfr 

3K 0 


(16) 


and Poisson’s ratio drops out. This equation is obviously correct when AT is just a 
constant; in this case AR = ( ay/3)ATRo , where ay/ 3 is the linear expansion coefficient. 

The equation (7) for AT(r ) contains the spherical Bessel function (6). The Bessel 
function is sin kr/kr, where k is complex and 

sin kr = sin ( qr - iqr) = cosh qr sin qr - /sinh qr cos qr . 

Here q = x/Rq. Using the above expressions in (16) results in the integrals 
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J r cosh qr sin q r ir - ^(sinh q r sin q r - cosh q r cos q r) + ^sinh v cos q r 


and 


f r sinh qr cos qr dr = — (cosh qr cos qr + sinh qr sin qr ) ^cosh qr sin qr , 

2 q 2 q 


ultimately giving 


AR _ a v (\-A v )R^F E 
4 K 


\ 2 =0 

— y.hj cos (jcp-Zj) 

“Plan ) H 


a v {l-A v )RlF E 
4 K 


CL \ 2 °° 

Y/*Ycos^cos£ ; + sin^sinQ . 

a Han ) h 


(17) 


In (17) 


hj=bj. 


Dr +D: 


D 


where 

D\ = [(1 + Co) x (sinh 2x - sin 2x)]/4x 2 , 

Di = {x (sinh 2x + sin 2x) - 2x 2 (cosh 2x + cos 2x) 

+ Co [cosh 2x - cos 2x - x (sinh 2x + sin 2x)] }/4x 2 , 

and 


(18) 


(19) 


( 20 ) 


Di = B 3 /2 . 


( 21 ) 
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Also, 


tan £; = 


A 

A 


(22) 


5. Ring particle along-track acceleration 


With the equation for AR in hand, the acceleration due to the thennal expansion 
and contraction can now be found. The acceleration of the particle due to solar radiation 
pressure is 


r 


x Rl C K F Sun 
cM s 


(23) 


where c is the speed of light, 

M= 4jrpR 3 /3 (24) 


is the mass of the particle, Cr is the radiation pressure coefficient, and F Sim is given by 
(1). For a diffusely reflecting sphere Cr » 1.2 (e.g., Smith et ah, 1991); this is the value 
adopted here. 

Assuming AR « Rq, so that R 2 ~ Ro 2 + 2Ro AR, the acceleration due to the change 
in radius is (dropping the subscript on R 0 ) 


$ 2 JtRARC R F Sun 

cM 


(25) 


where 


r s = cos s v x + sin e v z 
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is the unit vector in the direction from the planet’s center to the Sun, with ey being the 
solar elevation angle; Sy = 0 when the Sun is in the planet’s equatorial plane. 

It remains to find how thermal expansion affects the orbital evolution. For a 
particle in a circular orbit, the change in the semimajor axis a with time t is 


dci 2 .3 j y 

dt n 


(26) 


(e.g., Blanco and McCuskey, 1961), where 


C — f • r 

°EX v l EX 


is the along-track acceleration due to the expansion effect, with 


t = - sin <j) x + cos <)> y 


being the unit vector in the direction of motion. By the above and (25) 

^x= f Ex-t = ^^cos£ l/ 4R(F Sl „ I sin<))) , (27) 

CM 


where 


F Sun Sin <P = F E 





sin cp + 




sin 2(p 



1 sin 3 (j> + ... 


so that 
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F su„ sin (f> = ~r~ 


F E 

f a E ) 

2 J 

2 

\ a Plan j 

1 


2[( 1 + <5 lm)Vl-^ + l] sinW ^ 


m = 1 


(28) 


with d mj being the Kronecker delta: d mj = +1 when m = j and is zero otherwise. Plugging 
(17) and (28) into (27) yields terms with factors of the fonn sin jp sin imp, among others. 
Using 

sin jcp sin nup = [cos (j-in)cp - cos (j+m)(j)\/2 

and retaining only the first tenn on the right side for which j = m yields the secular terms 
in Sex, resulting in 


^ EX ~ 


3a v (] - A V )C R F 'l cos e 
32 Kpc 


\ a Plan / j - 1 


^ [(1 + Ay ) b j-i “ h j+i ] 11 j sin 


(29) 


The above expression gives the along-track acceleration averaged over one orbital 
revolution about the planet. It is interesting that the semimajor axis a does not appear in 
the coefficient multiplying the sum. It is implicit in the variables in the summation. 

The planet’s year must be averaged over next, because the tilt of the planet’s axis 
causes £y to constantly change. The bj in (29) are functions of £y and a; how much 
shadow is cast on the ring system depends on the solar elevation angle and the distance of 
the rings from the planet. The shadow cast on the equatorial plane is an ellipse, with a 
semimajor axis RpiaJsm £y and semiminor axis Rpi an , where Rpi an is the planet’s radius. A 
particle passes through the shadow when the circle of radius a and the shadow ellipse 
intersect. The intersection points can be found from the equations for the circle 



= 1 
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and the ellipse 


x 


2 


l Rpian \ 
\sin £ v ) 



= 1, 


yielding 


A = ■ 


2 Arc sin 

( R l,an- a2 S in 2 £ v )‘ /2 

a cos £ v 


L J 


2 n 


(30) 


which is to be used in (3) and (4). The solar elevation £y varies with time according to sin 
£v= sin £ sin (2nt/Ppi an ), where for example £ = 25° is the obliquity of Mars’ equator, and 
P plan = 1.6 yr is the length of the Martian year in Earth years. Due to the complexity of 
(30), at this point the analytical development is forsaken for a numerical one in order to 
average (29) over the planet’s year. 

The ring particle temperature To will be found from 


Ap crT = F 

^ t IR U1 0 r E 


( a E ) 

2 

\ a Pla„) 



Mi-M 




y PIan 

a 


£0-A„J(l -A m )(R 


Plan 


(31) 


The first term in the brackets is the average amount of sunlight shining directly on the 
particle. The second term gives the heating of the particle from visible light reflected off 
of the planet, where Api a „ the planet’s albedo in the visible (Conrath et ah, 1989); here it 
is assumed that the body reflects uniformly off its surface. The third term in the brackets 
gives the heating of the particle from infrared radiation coming up from the planet; Am is 
the particle’s albedo in the infrared. The parameter £ accounts for whether the planet 
emits more energy than it absorbs from the Sun. In the case of the terrestrial planets, it is 
assumed that % = 1, since the geothermal flux is much smaller that solar insolation. 
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However, f = 1.67 for Jupiter and 1.78 for Saturn. In other words, Saturn and Jupiter 
each have a significant source of internal energy (Conrath et ah, 1989). It is assumed in 
this tenn that the planet’s infrared radiation is emitted uniformly. 

It should be noted that To computed via (31) is not the same as the average 
temperature T avg obtained by measuring temperatures over space and time and taking the 
numerical average; instead T avg < To (Rubincam, 2004). The temperature dependence of 
K and C p is computed using To in the following sections. Using T avg would be a better 
choice, especially in the case of Saturn, where the ice particle temperatures may be close 
to 73 K, which is the temperature where the thermal expansion effect changes sign for 
polycrystalline ice (see section 9). However, the difference between To and T avg is ignored 
here, due to the complications of finding T avg . 

6. Ring particle orbital evolution 

The rings are assumed to be made of rock for the asteroids and planets, except for 
Saturn; there the rings are taken to be composed of water ice. For rock, each ring particle 
has a density of p = 3000 kg m' 3 and has a constant visible albedo Ay = 0.05, infrared 
albedo Am = 0.1, and infrared emissivity Em = 0.9. The volume coefficient of thermal 
expansion for rock is simply assumed to be a constant 

ay= 15 x 10" 6 

regardless of To, while the thermal conductivity K and specific heat C p are respectively 
assumed to depend on temperature the way basalt does: 

K = 3.0 - 0.002 To W m' 1 K' 1 

and 


C p = 500 + 0.833 To J kg' 1 K’ 1 
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(Roy et al., 1981, p. 434-435). Parameters for seven planets and two asteroids are given 
in Table 1. Further, in (29), the summation is cut off at j = 19. 

The Poynting-Robertson along-track acceleration Spr is also estimated here. It is 
taken to be 


where v = na is the speed of the particle in its circular orbit (Burns et al., 1979). The 


planet. Direct solar radiation pressure, the Yarkovsky effects, atmospheric and plasma 
drag, electric charging, magnetic fields, the solar wind, gravitational perturbations from 
other bodies (such as satellites), and the lumpy gravitation field of the body itself are all 
ignored. Only the Poynting-Robertson effect is considered here in addition to the thermal 
expansion effect. 

7, Rocky rings around Jupiter and the inner planets 

Figures 3-9 show da/dt for rocky particles in terms of planetary radii per 10 6 years 
as a function of distance from the orbited body, where dci/dt is the sum of the thermal 
expansion effect (29) plus the Poynting-Robertson effect (32): 


The figures show da/dt for four particle sizes: R = 0.001 m, 0.005 m, 0.01 m, and 0.02 
m. Sizes below 0.001 m are not shown because they give unstable results, apparently due 
to the number of terms used in the Fourier expansion of the shadow and to the 
approximate equations used here. Fewer Fourier terms increases the instability, as does 
being far from the orbited planet where the particles spend less time in shadow. 



(32) 


equation for Spr assumes the planet’s visible and infrared radiation is isotropic around the 


— = 2(S E x + S PR )ln . 


(33) 
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For the particle sizes and distances from the planet (1.1 - 7.0 Rpi an ) considered 
here, da/dt is positive for the inner planets; the thermal expansion effect is stronger than 
the Poynting-Robertson effect and the ring particles move outward. But for Jupiter and 
beyond the Poynting-Robertson effect can be stronger close to the planet for some 
particle sizes; in that case, da/dt < 0 and the rings move inward. 

Mercury (Fig. 3) shows the fastest ring particle orbital revolution in terms radii 
per million years compared to the other inner planets. The mean motion n in (33) does 
not show much variation for the inner planets, being roughly 0.001 N' 3 ' 2 s' for all of 
them, where N = a/Rpi an , so that the mean motion does not account for much of the 
differences in ring evolution in (33). Rather, Mercury is fast because its radius in small 
compared to Venus, Earth, and Mars, which appears to make it speedier in terms of the 
chosen units; and because it is closest to the Sun. It will be noted from (29) that the 
thennal expansion effect goes like square of the solar flux, i.e., FE 2 (a E /api a „) 4 , for the 
following reason: not only does the thermal expansion of a particle depend on the solar 
flux, but also the force on it also depends on solar radiation pressure (Fig. 1); hence the 
square. On the other hand, Poynting-Robertson depends on the flux coming up from the 
planet, which depends linearly on the light the planet receives from the Sun. Hence 
Poynting-Robertson varies only like the solar flux Fp(aE/api a „) 2 itself. For Mercury, the 
Poynting-Robertson effect, though also strong near the Sun, seldom accounts for more 
than ~15% of the thermal expansion effect for the sizes considered. 

Venus and Earth (Fig.s 4-5) have nearly the same radius and mass. Venus’s 
orbital evolution is -4 times the Earth’s, mainly owing to Venus being closer to the Sun. 
At the Earth, the Poynting-Robertson effect is -20% of the thermal expansion for 0.001 
m particles close to the Earth, and turns out to be larger than the thermal expansion effect 
at distances > 9 R Ear th (not shown in Fig. 5, which cuts off at 7 R Ear th )• 

Mars (Fig. 6) and Earth have nearly identical curves in terms of planetary radii 
per million years because while Mars is further from the Sun, which weakens the effect, 
Mars has only half the radius of the Earth, which appears to hasten its orbital evolution in 
the chosen units. The vertical dotted lines in Fig. 6 show the positions of Phobos and 
Deimos, which are possible sources of dust. At Jupiter (Fig. 7) the Sun is so weak that 
ring evolution is very slow, and the Poynting-Robertson effect with its high velocities 
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dominates close to the planet. The shaded areas in Fig. 7 show the approximate locations 
of the halo ring, the main rings, and the Amalthea gossamer ring. The Thebe gossamer 
ring is not shown. 

8. Rocky rings of Ida and an NEA 

The results for a hypothetical Near-Earth Asteroid (NEA) are shown in Fig. 8. 
The NEA is assumed to be a sphere with radius Rnea = 1000 m, density p = 2000 kg m" 3 , 
in the same orbit as the Earth (ci^ea = 1 AU), and with the same obliquity (23.5°). In this 
figure da/dt is measured in asteroid radii per thousand years, rather than per million 
years. The reason for the apparently quickened pace in Fig. 8 is the NEA’s small radius. 

Figure 9 shows the rates of orbital evolution for rocky rings around asteroid 243 
Ida. Ida is assumed to be a sphere here with a radius Rida = 1.57 x 10 4 m, but in reality the 
asteroid is highly elongated. Measured in terms of Ida radii per million years, da/dt 
appears fast due to the asteroid’s small size. 

9, Icy rings of Saturn 

Water ice can be polycrystalline, single crystals, or amorphous. Only the 
polycrystalline phase is considered here because the observed thermal inertia of the rings 
is so low. The ice is assumed to be “fluffy”, with a density p = 500 kg m' 3 rather than the 
918 kg m' 3 of a solid piece of ice. For polycrystalline ice K = 0.0001 W m' 1 K" 1 (Spilker 
et ah, 2003; Ferrari et ah, 2005), and C p = 715.0 + 8.2733 ( T 0 - 80.0) J kg' 1 K' 1 (Giauque 
and Stout, 1936). 

As for visible albedos used in (31), the observed albedos vary from ring to ring: 
for the C ring, 0.17 <Ay< 0.48, while for the B ring 0.41 <Ay< 0.73 and for the A ring 
0.31 <Ay< 0.63 (Porco et ah, 2005). Here the visible albedos Ay are simply taken to be 
0.35 in the C ring (1.1-1. 5 Rs at ) , 0.8 in the B ring (1. 5-2.0 Rs a t), and 0.5 in the A ring and 
beyond (2. 0-7.0 Rs a i), where Rs at is the radius of Saturn. The least-known parameters in 
(31) are the infrared albedo Am and infrared emissivity e m, these are simply set to Am = 
0.15 and £/r = 0.85 for all three rings. 
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The volume coefficient of thermal expansion ay can be found from Hobbs (1974, 
p. 347), who plots the observed linear thermal expansion of polycrystalline water ice as a 
function of temperature T (circles in Fig. 10). After converting Hobbs’ linear expansion to 
volume expansion by multiplying by a factor of 3, ay follows the approximate equation 

3 [0.0017 (To - 10.0) 2 - 6.75] x lO' 6 ^ 1 (34) 

for polycrystalline water ice. This equation is the solid line plotted in Fig. 10 and is used 
as the expression for ay. The line reproduces the polycrystalline curve quite well in the 
temperature range 30 - 130 K. The parameters were chosen so that (34) crosses ay = 0 
exactly at T 0 = 73 K, as does the observed polycrystalline curve. Hence (34) is a 
reasonable approximation for ay for the temperatures encountered at Saturn. 

When ay > 0 the particles are driven away from Saturn, as shown in Fig. 1. For 
temperatures below 73 K, Fig. 10 indicates that ay < 0 for polycrystalline ice, so that 
particles with low average temperatures would actually be driven toward Saturn and the 
rings would collapse. For sparsely populated rings the temperatures are near 73 K and the 
rings may expand or collapse from the thennal expansion effect, depending upon the 
exact temperature. 

Figure 1 1 shows the results for Saturn. In this figure To > 73 K everywhere. The 
approximate positions of the C, B, and A rings are denoted by the grey areas. The 
Poynting-Robertson effect dominates the thennal expansion effect for Saturn, and these 
rings would collapse from the combination of the two effects for the size range indicated. 

11. Meteoroids orbiting the Sun 

Meteoroids in eccentric orbits about the Sun will also suffer from the thermal 
expansion effect. The meteoroid’s radius will increase from solar heating as it approaches 
the Sun, but will not reach its maximum radius until after perihelion passage because of 
thermal inertia (Fig. 12). Far from the Sun the meteoroid cools and contracts, lessening 
solar radiation pressure. As in shadow passage, radiation pressure acts asymmetrically on 
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the orbit, increasing both semimajor axis a and the eccentricity e. The acceleration (25) 
will be of order 


H reaches a maximum of -0.0066 at radius -2 m for a rocky meteoroid in an orbit which 
takes it between Mercury and the Earth (details omitted). This leads to a secular rate of 


eccentricity. The eccentricity e itself will also change secularly, but the changes in both a 
and e are tiny. 

Dmitrievskii (1981) considers the effect thermal expansion has on the motion of a 
meteoroid’s perihelion, but assumes no thennal lags. Without the lags there are no 
secular effects on a and e. 

11. Discussion 

In computing orbital evolution, all rings are assumed to be sparsely populated: the 
ring particles are taken to be far enough apart so that mutual shadowing can be ignored. 
All ring particles are also assumed to stay in circular orbits in the equatorial plane of the 
body they orbit. Staying in such orbits for long periods of time would probably require 
the rings to be densely populated enough for mutual collisions to kill any developing 
orbital eccentricities or inclinations. Hence there is something of a contradiction: on the 
one hand the rings assumed to be sparsely populated, but on the other hand populated 
densely enough for collisions. 



by (17), where 



change in the semimajor axis of da/dt = - 0.0003 e 1 (AU/10 6 y), where e is the orbital 
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More serious is the issue of particle lifetimes. The other forces mentioned in 
section 6 which might influence the lifetimes are all ignored here. Hamilton and Bums 
(1992), for instance, find that solar radiation pressure is very efficient in increasing 
orbital eccentricities for small particles orbiting asteroids until the particles crash into the 
surface. Thermal stresses have been neglected. These stresses could lead to cracking and 
comminution of particles as they constantly plunge in and out of the planet’s shadow and 
their temperatures plummet and skyrocket. Any hysteresis in the expansion-contraction 
cycle is also ignored. Impacts, perhaps the most important destroyer of ring particles, 
have been neglected as well (Durisen, 1984; Durisen et al., 1996). 

The results indicate that the thermal expansion effect, though weak, outweighs the 
Poynting-Robertson effect for a limited range of sizes. As particles become smaller they 
become nearly isothennal and the thennal lags are tiny, so that the thennal expansion 
effect becomes negligible, while the rate of Poynting-Robertson drag increases. Thus for 
the smaller-sized particles, probably not much smaller than the lower limit of R = 0.001 
m investigated here, the Poynting-Robertson effect dominates. 

How does the sum of the thermal expansion effect plus Poynting-Robertson effect 
(TE + PR) compare with the sum of the seasonal Yarkovksy effect plus Yarkovsky- 
Schach effect (SY + YS)? While no detailed comparison is given here, for the terrestrial 
planets they are roughly equal at a = 2 Rp/ an when R = -0.02 m. At smaller sizes TE + PR 
dominates, while at larger sizes SY + YS dominates. (SY and YS for these planets were 
computed using the equations from Rubincam (2006); the seasonal effect was averaged 
over all spin axis positions.) Hence the thermal expansion dominates in the narrow niche 
-0.001 m to -0.02 m between Poynting-Robertson and Yarkovsky. 

Where might the thermal expansion effect be most significant, if at all? Ring 
particles have to exist for the thermal expansion effect to work on them. No rings are 
expected for Mercury and Venus because they lack an orbital source, and no rings have 
been observed. O’Keefe (1980) speculated that the Earth once had tektite rings; Hancock 
and Povenmire (2010) conjecture remnants may still survive. But the prospects for rings 
around the three innermost planets appear poor. 

NEAs might produce debris from Yarkovsky-0 ’Keefe -Radzieevski-Paddack 
(YORP) spin-up (Paddack, 1969; Rubincam, 2000), from YORP fissioning (e.g., 
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Jacobson and Scheeres, 2011), from tidal encounters with planets (e.g., Richardson et ah, 
1998), or from impacts. Moreover, the thermal expansion timescale is fast compared to 
the asteroid’s radius (Fig. 8), which should help orbital debris stay in orbit. Flowever, one 
suspects the timescales of solar radiation pressure (Hamilton and Bums, 1992) and the 
asteroid’s lumpy gravitational field are much faster still, overwhelming the thermal 
expansion effect for NEAs as well as for main-belt asteroids. Further, the Poynting- 
Robertson effect from an NEA is so weak that a fully relativistic treatment may produce 
solar Poynting-Robertson terms (Rubincam, submitted manuscript, 2013). 

Jupiter has rocky rings (e.g., Burns et al., 1984, 1999), but the planet is far from 
the Sun and its environment is not as benign for small particles as that of the inner 
planets. But it is intriguing that Poynting-Robertson drag dominates thermal expansion 
for the smaller particle sizes at the inner edge of Jupiter’s halo ring (Fig. 7). Inside the 
inner edge particles fall on Jupiter, while outside they are pushed away by the thermal 
expansion effect. Saturn’s icy rings have plenty of particles, but the planet is almost twice 
as far from the Sun as Jupiter, where the effect is already weak, so it is doubtful it plays 
much of a role at Saturn. 

Probably the best candidate is Mars, though so far no rings have been observed 
there (Showalter et al., 2006). Its moons Phobos and Deimos would supply the dust, and 
Mars is close enough to the Sun for the thermal expansion effect to be relatively strong 
(Fig. 6), helping keep rings alive. 

What about the thermal expansion effect operating on meteoroids orbiting the 
Sun? The effect requires an initial orbital eccentricity e to work. Due its weakness and the 
fact that it depends on e 2 , the effect is unlikely to have much application to meteoroids. 
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Table 1. Parameters for seven planets and two asteroids. 


Object 

@Plan 

(AU) 

Mpi a n 

(1 0 24 kg) 

Rpian 

(10 6 m) 

£ d 

(deg) 

Apian 


Mercury 

0.39 

0.330 

2.440 

1.0 

0.07 

1.0 

Venus 

0.72 

4.869 

6.052 

2.6 

0.90 

1.0 

Earth 

1.00 

5.974 

6.391 

23.5 

0.31 

1.0 

NEA 

1.00 

8.378 x 10 42 

0.001 

23.5 

0.05 

1.0 

Mars 

1.52 

0.642 

3.389 

25.0 

0.25 

1.0 

Ida 

2.86 

4.2 x 10' 8 

0.016 

1.0 

0.08 

1.0 

Jupiter 

5.20 

1898.6 

71.492 

3.13 

0.52 

1.67 

Saturn 

9.58 

568.8 

60.330 

26.73 

0.34 

1.78 


a Semimajor axis of the object’s orbit. 
b Mass. 
c Radius. 
d Obliquity. 
e Visible albedo. 

Ratio of energy emitted to that received from the Sun. 
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Figures 

Fig. 1. Schematic of the thermal expansion effect. A ring particle orbits around a body; 
the view is looking down on the particle’s orbital plane. Sunlight (wavy lines) causes 
radiation pressure on the particle (thick arrows). The size of the arrows indicates the 
magnitude of the force. Sunlight heats the particle, making it expand, increasing radiation 
pressure. Shadow passage (grey region) causes particle contraction, decreasing radiation 
pressure. The forces do not balance around the orbit because of thermal inertia, causing 
the orbit to expand secularly. 

Fig. 2. Schematic of a particle’s orbit. A particle P is in a circular, equatorial orbit 
around a body. Its distance from the center of the body is a. The Sun lies in the x-z plane, 
with elevation angle Sy. 

Fig. 3. Mercury ring particle rates of orbital evolution cla/clt as a function of distance for 
the combined thermal expansion effect and the Poynting-Robertson effect, for stony 
particles with radii of 0.001 m, 0.005 m, 0.010 m, and 0.02 m. The particles are assumed 
to be in circular orbits. The abscissa is distance from the center of the planet, measured 
Mercury radii R-Merc- The ordinate is given in terms of R.Merc per 10 6 years. The curves may 
differ slightly from their actual values due to the “quantum” nature of PowerPoint™ 
graphs. 

Fig. 4. Same as Fig. 3, but for Venus. 

Fig. 5. Same as Fig. 3, but for Earth. 

Fig. 6. Same as Fig. 3, but for Mars. 


Fig. 7. Same as Fig. 3, but for Jupiter. 
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Fig. 8. Graph similar to Fig. 3, but for a hypothetical spherical NEA with radius 1000 m 
and density 2000 kg m 3 and in Earth’s orbit with Earth’s obliquity. The ordinate is in 
radii per 10 3 years, instead of per 10 6 years as in the other figures. 

Fig. 9. Same as Fig. 3, but for asteroid 243 Ida, assuming Ida is a sphere. 

Fig. 10. The measured volume coefficient expansion for polycrystalline water ice as a 
function of temperature (circles) and the curve assumed here. The coefficient is negative 
for temperatures below 73 K, so that the thermal expansion effect would cause orbital 
decay rather than expansion for temperatures below this value. 

Fig. 11. Same as Fig. 3, but for ice particles orbiting around Saturn. 

Fig. 12 Schematic of the thermal expansion effect for a meteoroid in an eccentric orbit 
about the Sun. The size of the arrows indicates the magnitude of solar radiation pressure 
on the meteoroid. As it approaches the Sun, sunlight (wavy arrows) heats up the 
meteoroid more and more, making it expand, increasing its cross-section and thus the 
force. The meteoroid achieves maximum size after perihelion passage, thanks to thermal 
inertia. Far from the Sun the weakened sunlight causes the meteoroid to contract, 
decreasing radiation pressure. The forces do not balance around the orbit, causing the 
orbit to expand secularly. The thermal expansion effect will also secularly change the 
orbital eccentricity. 
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Fig. 1 . Schematic of the thermal expansion effect. A ring particle orbits around a body; 
the view is looking down on the particle’s orbital plane. Sunlight (wavy lines) causes 
radiation pressure on the particle (thick arrows). The size of the arrows indicates the 
magnitude of the force. Sunlight heats the particle, making it expand, increasing radiation 
pressure. Shadow passage (grey region) causes particle contraction, decreasing radiation 
pressure. The forces do not balance around the orbit because of thermal inertia, causing 
the orbit to expand secularly. 
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Fig. 2. Schematic of a particle’s orbit. A particle P is in a circular, equatorial orbit 
around a body. Its distance from the center of the body is a. The Sun lies in the x-z plane, 
with elevation angle £y. 
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0.001 



Fig. 3. Mercury ring particle rates of orbital evolution da/dt as a function of distance for 
the combined thermal expansion effect and the Poynting-Robertson effect, for stoney 
particles with radii of 0.001 m, 0.005 m, 0.010 m, and 0.02 m. The particles are assumed 
to be in circular orbits. The abscissa is distance from the center of the planet, measured 
Mercury radii Rmck- The ordinate is given in terms of RMerc per 10 6 years. The curves may 
differ slightly from their actual values due to the “quantum” nature of PowerPoint™ 
graphs. 
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Fig. 4. Same as Fig. 3, but for Venus. 



Fig. 5. Same as Fig. 3, but for Earth. 
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Fig. 6. Same as Fig. 3, but for Mars. 
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Fig. 7. Same as Fig. 3, but for Jupiter. 
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Fig. 8. Graph similar to Fig. 3, but for a hypothetical spherical NEA with radius 1000 m 
and density 2000 kg m' 3 and in Earth’s orbit with Earth’s obliquity. The ordinate is in 
radii per 10 3 years, instead of per 10 6 years as in the other figures. 



Fig. 9. Same as Fig. 3, but for asteroid 243 Ida, assuming Ida is a sphere. 
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Fig. 10. The measured volume coefficient expansion for polycrystalline water ice as a 
function of temperature (circles) and the curve assumed here. The coefficient is negative 
for temperatures below 73 K, so that the thermal expansion effect would cause orbital 
decay rather than expansion for temperatures below this value. 



Fig. 11. Same as Fig. 3, but for ice particles orbiting around Saturn. 
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Fig. 12 Schematic of the thermal expansion effect for a meteoroid in an eccentric orbit 
about the Sun. The size of the arrows indicates the magnitude of solar radiation pressure 
on the meteoroid. As it approaches the Sun, sunlight (wavy arrows) heats up the 
meteoroid more and more, making it expand, increasing its cross-section and thus the 
force. The meteoroid achieves maximum size after perihelion passage, thanks to thermal 
inertia. Far from the Sun the weakened sunlight causes the meteoroid to contract, 
decreasing radiation pressure. The forces do not balance around the orbit, causing the 
orbit to expand secularly. The thermal expansion effect will also secularly change the 
orbital eccentricity. 



